Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Pulsed radiation from neutron star winds 
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Abstract. The radiation of a pulsar wind is computed assuming that at roughly 10 to 100 light cylinder radii 
from the star, magnetic energy is dissipated into particle energy. The synchrotron emission of heated particles 
appears periodic, with, in general, both a pulse and an interpulse. The predicted spacing agrees well with the Crab 
and Vela pulse profiles. Using parameters appropriate for the Crab pulsar (magnetisation parameter at the light 
cylinder (tl = 6 x 10*, Lorentz factor F — 250) agreement is found with the observed total pulsed luminosity. This 
suggests that the high-energy pulses from young pulsars originate not in the corotating magnetosphere within the 
light cylinder (as in all other models) but from the radially directed wind well outside it. 
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1. Introduction 

Recently, the non-axisymmetric "striped pulsar wind" in- 
vestigated by Coroniti (199C|) and Michel (1994) has been 
reexamined (Kirk & Lyubarsky, 2001). A striped wind 
is a structure similar to a Parker spiral, produced by 
the outward radial advection of magnetic field lines with 
foot points anchored on the surface of the rotating neu- 
tron star. If the field lines originate from the poles of an 
obliquely rotating dipole, there exists a region around the 
equatorial plane where the polarity of the field at a fixed 
radius reverses at the rotation period, producing stripes 
of alternating magnetic field direction. Magnetic energy 
is dissipated into parti cle energy if rec onnection occurs at 
the stripe boundaries. [Coroniti (199C ) proposed that this 
kind of dissipation proceeds at a rate sufficient to main- 
tain the thickness of the current sheet appro ximately equa l 
to the gyro-radius of the heated particles. Michel (1994 ) 
used essentially the same criterion, phrased in terms of 
the velocity of the current carriers. These early papers 
concluded that observations of the Crab Nebula, which 
indicate conversion of Poynting flux to particle-borne en- 
ergy flux within IO^tl, (where tl is the radius of the light 
cylinder) could be explained by this mechanism. However, 
Lyubarsky & Kirk (2001) point out that reconnection is 
accompanied by acceleration of the wind. Both processes 
occur on the same timescale in the Coroniti/Michel ap- 
proach, which is also the timescale on which the plasma 
expands. As a result, time dilation in the accelerated flow 
significantly reduces the dissipation rate. In the case of the 
Crab, conversion of Poynting flux to particle-borne energy 



flux is not achieved until IO^^tl — to o slow to accou nt for 
the observations [the cr-problem; see Melatos (1998| )]. 



However, our present understanding of the physics of 
the dissipation process does not exclude the possibility 
that it proceeds much more rapidly than suggested by 
Coroniti and Michel. Recent investigations of reconnection 
in a relativistic pair plasma indicate that the timescale is 
related to the crossing time of a fast magnetosonic wave 



(Zcnitani & Hoshino, 2001). In the case of a transsonic 
wind, the process could then be completed within a few 
pulsar periods, giving the plasma insufficient time to ac- 
celerate to high Lorentz factor. If rapid reconnection in 
fact occurs, it will be accompanied by radiation losses 



(Melatos, 1998) that may give rise to an observable signa- 
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ture. In this letter we adopt this hypothesis. We assume 
that dissipation is triggered at a surface r — tq well out- 
side the light cylinder in the supersonic MHD wind of a 
pulsar and takes place rapidly. Using two simple models 
of the electron distribution we show (i) that for the winds 
of young pulsars (e.g., the Crab) the predicted emission 
is pulsed, with, in general, a pulse and interpulse that 
are not symmetrically spaced in phase and (ii) that the 
predicted synchrotron luminosity in the pulses agrees to 
order of magnitude with that observed in the high energy 
emission of the Crab pulsar. 

Our calculations suggest that the high energy pulses 
emitted by young pulsars originate as synchrotron emis- 
sion in the pulsar wind, far outside the light cylinder. This 
is in marked contrast with conventional models of high- 
energy pul se production at an inner or outer gap [for a 
review see Harding (2001 )] and also with magnetospheric 



synchrotron based models ( Machabeli et al, 200C ; Crusius 



2 



J. G. Kirk et al.: Pulsed radiation from neutron star winds 




Fig. 1. A meridional section of the current sheet of an 
obHque rotator with a radial wind. The thin lines limiting 
the extent of the sheet in latitude depict the plane of the 
magnetic equator when the magnetic axis lies in the plane 
of the section. 



Watzel et al, 2001). The current sheet just beyond the 
region of closed field lines in the magnetosphere has been 
proposed as the source of high energy radiation from the 
Crab pulsar (Lyubarsky, 1996). But, in common with the 
others, this model relies for pulse production on the coro- 
tation with the star of a, bea.m of radiation emitted frmn 



(i) The emission region is located on the current sheet in 
the supersonic MHD wind of a rotating neutron star 
with an oblique, split-monopole magnetic field. 

(ii) The dependence on frequency v of the volume emissiv- 
ity is a power law: e oc 

(iii) The emission switches on abruptly as plasma in the 
sheet crosses the surface r = rg . Thereafter the volume 
emissivity decays as a power law in radius: e oc r~^~'. 

Outside the light cylinder, assumption (i) should be a 
good approximation, independent of the actual configura- 
tion of the magnetic field close to the pulsar surface. The 
assumption that emission occurs only at the exact posi- 
tion of the current sheet is an idealisation — in reality, the 
entire flow may be modulated, making the observed pulse 
more complex. In this case, our computations provide the 
Green's function with which the pattern is to be folded to 
find the observed pulse shape. If the wind carries enough 
particles to ensure MHD behaviour at least i nitially [for an 
alternative approach see Michel fc Li (1999| )] , a simple so- 
lution is available in the ultra-relativistic limit (Bogovalov, 



within Ithe light cylinder. 



2. Geometry of the emission region 

Any radial, relativistic flow containing a periodic mod- 
ulation of the emissivity is likely to appear pulsed to an 
observer because of two effects. Firstly, the strong Doppler 
boosting of approaching parts of the flow means that only 
a small cone of the flow is visible, propagating within an 
angle of roughly l/F with respect to the line of sight, 
where F = (1 — /3^)~^/^ is the Lorentz factor and c/? the 
radial speed. The spread in arrival times at the observer 
of photons emitted at the same radius but from different 
parts within this cone is roughly r/ (cF^). Secondly, if emis- 
sion is restricted to a range Ar in radius, the fact that the 
photons move only slightly faster than the flow when seen 
in the laboratory frame means that the spread in arrival 
times is roughly Ar/(cF^). Thus, for a wind modulated at 
the period P = 27rrL/c of a rotating star, the approximate 
condition for the observation of pulsed emission is 



Ar, r ^ F^tl 



(1) 



If this is fulfilled, the observed pulse shapes appear simi- 
lar to the pattern of modulation of the flow. In particular, 
for a flow containing the corrugated spiral current sheet 
of an oblique rotator (Fig. Q), there should in general be 
two components (pulse and interpulse) which are symmet- 
rically placed if the pulsar is viewed from the equatorial 
plane, but are asymmetrical when viewed from higher lat- 
itudes. Denoting by a the angle between the rotation axis 
and the magnetic axis, and by C the angle between the ro- 
tation axis and the viewing direction, these pulses appear 
only for |90° - a| < C < 180° - |90° - a\. 

To improve on the estimate (0) and calculate the ex- 
pected pulse shapes we adopt the following model of the 
volume emissivity e of the wind: 



1999). In it, the velocity is radial and constant and the 



magnetic field is almost toroidal. The thin current sheet 
in such a wind is described by the equation 



r-s(e,0, t) 

Pri^ [± arccos (— cot a cot ff) + ct/ri^ 



2mr] (2) 



( Kirk fc Lyubarsky, 2001 ), where {r^9,4>) are polar coor- 
dinates, t is the observer-frame time and n is an integer. 

Assumption (ii) ensures the pulse shape is indepen- 
dent of frequency. In reality, the electron distribution may 
exhibit breaks and cut-offs, which translates into differ- 
ent pulse shapes at different frequencies. Our assumption, 
however, ensures that all parts of the flow can contribute 
to the emission at a given frequency, making it harder to 
obtain sharply defined pulses. 

The dependence of the observed flux F on t, a and the 
direction of observation n is: 

/+00 />oo />90"+a 

dt' dr d9 
-oo Jro J90°-a JQ 

^-q^-i+a^ [r - n{6, 0, t')] S{t' -t + n-r/c) (3) 

where D — 1/ (1 — /?r ■ n) is the conventionally defined 
Doppler factor, apart from a (constant) factor of F. 
Choosing a = and q = 3, which is appropriate for the 
synchrotron radiation of a fixed number of particles that 
undergo predominantly adiabatic losses in the expanding 
plasma, we present sample pulses for various Lorentz fac- 
tors in Fig. 1^. For the higher Lorentz factors this figure 
also displays the effect of a faster fall-off of the emissivity 
towards higher radii {q — 6), mimicking the effect of an 
accelerating wind. The asymmetry in time of the pulses 
is produced by our assumption of a sudden switching on 
of the emission and is reduced by a more rapid switch- 
ing off. The parameters of this figure have been chosen to 
correspond to the Crab pulsar. The radio emission of this 
object is complex. The shape of the high energy pulses 
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Fig. 2. Pulse shapes obtained by a numerical integration 
of Eq. (^. Four values of the Lorentz factor (F) are shown, 
the remaining parameters are: a = 60°, C = 60°, tq = 
30rL, a = and q — 3 (see text). The total emission is 
shown as a solid line. For F = 2 and 5, the contributions of 
the current sheets of different polarity are shown as dotted 
lines. For F = 10 and 50, the dotted line shows the effect 
of a more rapid fading of the emissivity to larger radius 
(9 = 6). 



is repeated in the radio, suggesting a common emission 
site, which in our model would be the wind. But there ar e 
additional radio components ( Moffatt fc Hankins, 199^) 



one of which displays 'core emission' properties ( Rankin, 
1990|)~|This suggests emission from close to the polar cap, 
with the magnetic axis passing within a few degrees of 
the line of sight. The core emission is in the precursor 
of the main pulse, which constrains the relative phase of 
the polar cap emission and wind emission in this inter- 
pretation. However, a complete solution linking the inner 
magnetosphere and the wind structure would be required 
to check this point. Nevertheless, core emission indicates 
that the inclination angle a of the magnetic axis to the 
rotation axis roughly equals the angle C between the rota- 
tion axis and the line of sight. The latter is determined by 
X-ray and optical observations of the 'torus' (Aschenbach 
& Brin |kmann, 1975| ; [Hester et al, 1995D to be 60°. In Fig. ^ 
the emission has been assumed to switch on at vq = 30rL 
and pulses appear for F > 5, in rough agreement with 
the estimate (|^). At higher Lorentz factors, the pulses 
sharpen. From Eq. (|^) it can be seen that the emission 
splits naturally into two contributions, corresponding to 
different signs of the term arccos(— cot a cot 9). In the first 
two plots of Fig. H these contributions are plotted sepa- 
rately as dotted lines. Although the detailed pulse shapes 



will depend on the structure of the reconnecting sheets in 
the Crab wind, as well as the details of the switching on 
mechanism and the amount of acceleration of the wind, 
the separation of the two peaks is not sensitive to these 
unknowns, being determined by the location of the sheets. 
This separation in phase equals arccos(cotacot(^)/180°, 
corresponding to 0.39 for the Crab pulsar, in good agree- 
ment with observation [e.g., Kanbach (1998); Kuiper et a] 
"(2001D ] 



For the Vela pulsar, conflicting interpretations of the 
polarisation sweep and the X-ray morphology give a « 
71°, C ~ 65° dRadhakrishnan fc Deshpande, 2001] ) or a w 
C Ki 53° (Helfand et al, 2001). These imply separations of 
the gamma-ray (sub-)pulses of 0.45 and 0.31 of a period 
respectively, so that the observed spacing of 0.4 marginally 
favours the former interpretation. 

3. Luminosity estimates 

Electrons heated in the current sheets will enter the sur- 
rounding magnetic field and emit synchrotron radiation. 
An estimate of the power radiated can be found by as- 
suming a monoenergetic electron population of Lorentz 
factor 7 (^ 1) and requiring pressure equilibrium between 
the hot plasma and the surrounding magnetic field. In the 
delta-function approximation for the synchrotron emissiv- 
ity, the luminosity at frequency into a solid angle f2 is 



dL/dt/dl7 « / dr AVn'^Ps.pS {ly 

J rn 



(4) 



(Skjaeraasen & Kirk, 2001), where A (< 1) is the fraction 



of a wavelength of the wind pattern occupied by hot elec- 
trons of proper density n'^, Pg.p. = (4/3)c(TT7^(i?'^/87r) is 
the synchrotron power emitted by a single electron, and 
z/Q is the characteristic frequency of emission, measured 
in the observer frame (ctt is the Thomson cross section). 
Pressure equilibrium implies n'^^mc^ = i3'^/87r, where 
B' is the magnetic field in the comoving frame and m the 
electron mass. Well outside the light cylinder, but before 
energy release begins, i.e., at tl < r < tq, the density 
and magnetic field scale as and respectively, and 
can be characterised by the values riL and Bl obtained by 
extrapolating back to r = tl- Provided the flow does not 
accelerate significantly whilst emitting synchrotron radi- 
ation, these same scalings, which imply 7 independent of 
r, can be used in the integration in Eq. (^). This is con- 
venient for an estimate of the luminosity, but is sensitive 
to the assumption that the wind docs not accelerate upon 
reconnection. 

If we identify the ratio of Poynting flux to kinetic 
energy (at the light cylinder) ctl = B^/{AnTni,mc^), 
define an average fraction of hot particles A = 
^0 X-Q drA/r^ ^ 1, and use the standard relation i?L = 
2.9 X io8pV2p-5/2 ^jjj, luminosity becomes 



dL 



= 3.9 X 10' 



'crArL-P\ 

FVoP* ] 



{P in sees) 



(5) 
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where dL /dQ is the synchrotron himinosity per steradian 
divided by the wind luminosity per steradian. Similarly, 
the characteristic frequency is 

hvo{r) = 2.6 X 10"^ (rL/r) cr^F^/^P^^/^ MeV. (6) 

The spectrum of this model with monoenergetic elec- 
trons extends to a maximum frequency given by Eq. 
with r = ro and is flat: dL/di/d^l oc z/". However, this 
aspect is sensitive to expansion and acceleration of the 
plasma after reconnection. If, in a more realistic scenario, 
the electron density is a power law iV(7) oc with 
p > 2, then most of the energy resides in the particles 
with lowest 7 and the above estimate (^) refers to the 
photons emitted by these particles. The emission at fre- 
quencies < i'o(''o) is dominated by the same low energy 
particles radiating at larger radius in a weaker magnetic 
field, and is determined by the expansion and accelera- 
tion of the hot plasma. In the absence of acceleration, 
the fiat spectrum derived above prevails. However, the in- 
trinsic spectrum, which is a power law dL/dvdQ oc i>~°', 
with a > 0.5, is revealed at frequencies v > i'o('"o): where 
particles radiating at tq dominate. In the case of the 
Crab pulsar, the optical to X-ray emission is consistent 



with a flat spectrum (Shearer & Golden, 2001), which 
steepens to a > 1 at a photon energy of approximately 
1 MeV ( iKanbach, 199^ ; [Kuiper et al, 200l| ). Equation (|[) 
then implies a = 1.1 x W^^/ro/ri^. Observations suggest 
dL/dil w 10~^, which, from Eq. (|^), is consistent with a 
mildly supersonic flow, F « 250, that starts to radiate at 
ro = 30rL. This satisfies the lower limit of 5rL (Bogovalov 



& Aha:onian, 2000) on the radius at which Poynting fiux 
can be converted to particle-borne energy flux. 

4. Discussion 

Current models of the high energy emission from rota- 
tion powered pulsars fal l into two group s: "polar cap" and 
"outer gap" models [see Harding (2001 )]. Each locates the 
emission region in the corotating magnetosphere of the 
star. The predicted pulse shapes depend sensitively on 
the uncertain geometry of the magnetic field in this re- 
gion, usually assumed dipolar. In contrast, in our model, 
pulses are emitted outside the corotation region. If the 
pulsar drives a supersonic, MHD wind, the field geometry 



there approaches a simple asymptotic solution ( Bogovalov, 



1999), determining the basic properties of the pulses. 

Pacini {Wff) and [Pacini fc Salvati (1983| , |1987| ) sug- 
gested that particles close to the light cylinder emit syn- 
chrotron radiation and are responsible for the optical 
pulses of the Crab and other pulsars, a theory which 
appears to be in reasonable agreement with observation 
( ^hearer fc Golden, 200l| ). Our model is similar in that 
the radiation mechanism is synchrotron emission and the 
scalin g of Eq. (|^) is close to that originally given by Pacini 



1971 ). It goes further, however, by specifying for the emis- 
sion region a precise geometrical structure and location. 
It is known that the dilution of the plasma in the 



wind must lead to non- ideal MHD behaviour ( Usov, 1975 ; 



Michel, 1982) and so can be responsible for triggering the 
emission. The position at which this occurs depends on the 
initial concentration of charges in the current sheet and is 
taken as a free parameter. The major uncertainty in the 
model is the speed with which dissipation proceeds. Our 
computations implicitly assume rapid dissipation (over a 
scale small compared to the radius). In this case, it is 
possible to make a rough estimate of the synchrotron lu- 
minosity which is in agreement with observations. It is 
conceivable that at least part of the radio emission could 
be produced in the wind region, as is suggested by the sim- 
ilarity of the pulse profiles at all frequencies seen in the 
Crab pulsar. However, until a candidate coherent mecha- 
nism can be identified, this connection remains very spec- 
ulative. Nevertheless, the wind scenario is a viable alter- 
native to current theories of gamma-ray. X-ray and optical 
pulses from rotation-driven pulsars and is capable of mak- 
ing testable predictions about the pulse profile. 
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